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We consistently analyze two-pion correlation functions and three-pion correlation functions of 
130A GeV Au+Au collisions measured at RHIC, applying three models of partially coherent pion 
sources. The effect of long-lived resonances on the chaoticity in the two-pion correlation function is 
estimated on the basis of a statistical model. We find that the chaoticity extracted from the three- 
pion correlation function is consistent with that extracted from the two-pion correlation function 
at vanishing relative momenta after the subtraction of the contribution of the long-lived resonance 
decays. All three models indicate that pions are emitted from the mixture of a chaotic source and 
coherent sources in the 130^4 GeV Au+Au collisions. 

PACS numbers: 25.75.Gz 



I. INTRODUCTION 



Pion interferometry is one of the most powerful tools in relativistic heavy ion physics. The fact that two-boson 
intensity correlation functions provide information concerning the geometry of particle sources is known as the Hanbury 
Brown-Twiss (HBT) effect. The two-pion correlation function C 2 (pi,p 2 ) has been analyzed extensively to explore 
space-time structures of collision dynamics pj. The HBT effect is based on the assumption of a chaotic source; i.e., 
particles are emitted with relatively random phases. The chaoticity index A in the two-particle correlation function 
is defined as the two-pion correlation strength, A = C^p, p) — 1. It is unity for a chaotic source and vanishes for a 
coherent one. Thermal emission is naturally regarded as a perfect chaotic source, @>I3 but non-thermal components 
need not be perfectly chaotic. For example, if a disoriented chiral condensation (DCC) region is produced, this 
domain can decay into coherent pions. Therefore, the chaoticity index A is an important parameter for understanding 
multiple particle production. However, experimental data are not so simply understood; even for a chaotic source, 
the chaoticity can be reduced by some other effects, such as long-lived resonances 0, IS El 01 HI ■ 

As an alternative tool to investigate particle sources, the three-particle correlation function has also been proposed 
H IIS IH HI Regarding the chaoticity of the sources, the three-particle correlation is a more promising tool than the 
two-particle correlation, because long-lived resonances do not affect a normalized three-pion correlator (third order 
cumulant correlation function), 
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Here, C 2 (pi,p 2 ) and C 3 (pi,p2,Ps) are the two- and three-pion correlation functions defined as 

C 2 {Pi,P2) = tj77 — ttttt — V' 2 

Wl{p 1 )Wl{p2) 

and 

n i \ W 3 (p 1 ,p2,p 3 ) , Q , 

C 3 Pl,P2,P3 = TT7-, 77777 V' 3 ) 

W 1 {p 1 )W 1 (p 2 )W 1 (p 3 ) 

with W n (p±, ■ ■ ■ ,p n ) being the n particle distribution. The chaoticity index in the three-pion correlator is the weight 
factor oj — r${p,p,p)/2, which becomes unity for a chaotic source. The three-pion correlator used in the current 
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experiments is a modified form of Eq. (JTJ, introduced due to a lack of statistics, as 

[C 3 (Q 3 ) - 1] - [C 2 (Q 12 ) - 1] - [C 2 (Q 23 ) - 1] - [CaCQsi) - 1] 



r 3 (Q 3 



- 1][C 2 (Q 23 ) - l][C 2 (Q 3 i) - 1] 



(4) 



Q3D and w is denned as ra(0)/2. In both definitions, 



where Q tj = yj-{pi ~ Pj) 2 and Q 3 = \/Q\ 2 + Q23 
expected to approach the same value at small Q 3 . 

In this paper, we analyze both the three- and two-pion correlation data measured by the STAR collaboration |l3j 
for central events in Au+Au collisions at 13QA GeV using three kinds of partially coherent models. The STAR 
collaboration, adopting a partial coherent model for the pion source, concludes that about 80% of pions are emitted 
from a chaotic source in the central collision events, and the chaoticity depends on the centrality In Ref. [lflj . 
analysis is carried out only with the three-pion correlation function, because the two-pion correlation function suffers 
from the effect of long-lived resonance decays. Hence, the chaoticity extracted from the three-pion correlation differs 
from that measured as the strength of the two-pion correlation function. One of the aims of this work is to account 
for this discrepancy. After an appropriate correction for the long-lived resonance decay contribution to the two-pion 
correlation function, the "true" chaoticity of the two-pion correlation function thereby obtained should be consistent 
with that extracted from the three-particle correlation function. We show that the chaotic fraction and the number 
of coherent sources evaluated from the "true" chaoticity of the two-pion correlation function are consistent with those 
obtained from the three-pion correlation function. Then, we study the structure of the source using three kinds of 
partial coherent source models. This paper is organized as follows. We give a brief explanation of the three kinds 
of source models in the next section. In Sec. IIIII we make corrections to the two-pion correlation for the long-lived 
resonance decay contributions. Section IIVI is devoted to a combined analysis of two- and three-pion correlation 
functions and analyses using the three models. 




FIG. 1: Schematic depiction of Model I-III. The shaded areas represent the chaotic sources, 
represent the coherent sources, (a) Model I, (b) Model II, (c) Model III. 
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II. MODELS OF PION EMISSION 



Let us start by giving a brief explanation of the models. In Model I, a chaotic source and a coherent source are 
assumed to be produced in every event of the collisions. This is a very limited case, because, in general, more than one 
coherent source can be produced in a collision. In Model II, we consider the case in which multiple coherent sources 
are produced. Model III is a mixture of Models I and II. A schematic depiction of the models is given in Fig. 

Model I is a partially coherent model |l(J in which the solitary parameter e pc represents the ratio of the particle 
number emitted from the chaotic source to that emitted from both sources. In this model, pions are emitted from a 
mixture of a chaotic (thermal) and a coherent source. Here we do not discuss the origin of the coherent component. 
In general, we need to fix the source function of the coherent component, as well as that of the chaotic one, in order 
to express the correlation functions. Fortunately, because the chaoticity A and the weight factor u> are given by the 
correlation functions at vanishing relative momenta, they can be expressed in terms of the chaotic fraction as 



A = e pc (2 



-pcj 



(5) 



Model II is a multicoherent source model, |l2j | in which the source is composed of a large number of similar, inde- 
pendent coherent sources. Here we do not discuss the origin of the coherence. Because each small coherent source is 
independent, a chaotic source is realized as a cluster of an infinite number of the coherent sources. Here, A and ui are 
given by 



to 
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where a m represents the mean number of coherent sources, which obeys the Poisson distribution. Note that each of 
the above two models contains only a single parameter (e pc for Model I and a m for Model II) for two quantities, the 
chaoticity A and the weight factor u>. Model III is a mixture of these two [see also Fig.^c)] and has already been 
applied to the SPS data by one of the present authors (H. N.). This model contains two parameters, which are 
related to A and u> as 

a + (1 — e) z 

_ 2a 2 + 2a(l-e) 2 + 3(l-£) 3 (l-2e) / a + (1 - e) 2 
U ~ 2[a 2 + 3a(l-e) 2 + (l-e) 3 ] V a ' () 

where e and a are the ratio of the number of particles from the chaotic source to the total number and the mean 
number of coherent sources, respectively. 

Substituting A and lu from the experimental data into the above equations, we can obtain the parameters of the 
models, e and a. In Models I and II, A and lj must provide consistent values for the single parameter e pc or a m . 
As noted in the previous section, the STAR group analysis of the weight factor uj in the central Au+Au collisions at 
y/SNN = 130 GeV based on Model I indicates a value of e pc — 0.8, which corresponds to A ~ 0.94 [13j . Such a large 
A seems to contradict the two-particle correlation data. However, we make a correction for the long-lived resonances 
in the two-pion correlation function before proceeding to conclusive calculations to obtain the parameters e and a 
using Eqs. (|5j)-(|3J|. Then we show that the analysis of STAR is consistent with the value of e pc obtained from the 
resonance-corrected chaoticity. 

III. LONG-LIVED RESONANCE DECAY CONTRIBUTIONS TO A 

The existence of long-lived resonances that decay into pions causes the apparent reduction of the chaoticity in the 
two-pion correlation function. This is because the pion source caused by the long-lived resonance decay spreads over 
a much larger scale than that which is measurable with the present experimental resolution. The reduction factor A cff 
is related to the fraction of pions from long-lived resonances as 

with Nn being the total number of emitted pions and being the number of pions emitted from long-lived reso- 
nances. We estimate the fraction of the tt~ originating from the resonances with the help of a statistical model. 
According to the recent analyses, this statistical model provides a simple and good description of the particle 
yields of the RHIC experiments (l5L Il6| . We take account of the resonances up to S*(1385), in accordance with 
Ref. @- We treat K®, r/, 77', 0, A, E and 5 as the long-lived resonances whose widths are smaller than 5 MeV. Heav- 
ier resonances can be ignored, because their contributions are thermally suppressed. The thermodynamic param- 
eters are determined with a ^ 2 fitting to the experimental data for the particle ratio. Here we take into account 
p/p,p/ir~ ,K~ /n~ ,K*°/h~ ,K*°/K*°, A/A and H/H. From the fit, we obtain the temperature T = 158±9 MeV and 
the chemical potential with respect to the baryon number /ib = 36± 6 MeV, with x 2 /dof = 2.4/5 ^tJ- For simplicity, 
we fix the chemical potential with respect to the third component of the isospin as fj,i 3 = 0. This does not affect 
the result, because the experimental data show ir~ /ir + ~ 1. The ratio of the number of particles i to the number of 
particles j can be calculated using the number densities in the local rest frame 0] . 

Then, the true chaoticity A truc , which is characterized by the partial coherence of the source, is estimated as 
A true = A exp /A eff , with use of Eq. 10. In general, A cxp obtained from the 1-dimensional Gaussian fitting can differ from 
that extracted from the 3-dimensional Gaussian fitting because of the projection, experimental resolution, and other 
effects, though these should be the same for an ideal measurement. In this paper, we use the value extracted from 
the 3-dimensional Gaussian fitting 

C 2 fit (q) = 1 + A cxp exp (-i? 2 idc <z s 2 idc - R 2 out q 2 out - K ng qL g ) , (10) 

as A Gxp . The STAR collaboration measured 0] A cxp for three transverse momentum bins. Because the transverse 
momentum dependence of A exp agrees with the results obtained from Eq. © Q , A cxp for the entire momentum range 
can be obtained by computing an average of A^ xp in which the square of the number of pions is used as the weight: 

r x P = EUK^f^un^dhidN/hdhr 
ELiL th un k tdh(dN/k t dk t r 
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Here, A^ xp and dN/k t dk t denote the measured chaoticity in the i-th bin and the pion transverse momentum distribu- 
tions, respectively. We give a derivation of Eq. in the appendix. In the calculation of Eq. (|llfl . we use a result 
of the Bose-Einstein fit to the pion spectra carried out by the STAR collaboration [23] • 
Using Eq. dHJ , we find A° XP = 0.57 ± 0.06. 

Combining A cff obtained from the statistical model calculation, we finally obtain A true = 0.93 ± 0.08 for the present 
parameter region. The uncertainty of 0.08 is the sum of the experimental uncertainty on A° xp calculated from 
uncertainties on each A^ xp and dN/k t dk t and the statistical error propagated from the fit of the thermodynamic 
quantities T and /j,b at 1-a confidence level. 



IV. COMBINED ANALYSIS OF THE TWO- AND THREE-PARTICLE CORRELATION FUNCTION 

In order to obtain the weight factor ui, we have to extrapolate r 3 (Q 3 ) [Eq. (QJ] to Q% = 0. Because r 3 (Q 3 ) is 
constructed from C2(Qij) — 1 and Cz(Qi) — 1 [Eq. |@J|], we need to establish the functional forms of C^iQij) and 
Cz{Qz) so as to reproduce the experimental data (Figs. Eland EJ- Although quadratic and quartic fits to r 3 (Q 3 ) at 
low Q 3 are used in Ref. [l3j, we adopt an alternative approach, because the power law fits are relevant only for small 
Qs, and it should be better to extrapolate ^(Q^) to Q3 = after the fits of C2 and C3 within finite but broad relative 
momentum ranges. We also note that the Q3 dependence of r 3 may be more complicated, [TlT ] and the quadratic and 
quartic fitting tends to lead to an overestimated weight factor. In the following, we treat a set of n~ data measured in 
the lowest momentum bin for the event with the highest multiplicity. We construct C2 and C3 from simple source 
functions with a set of parameters. This construction is determined by minimizing x 2 with respect to the experimental 
data. For simplicity, we used a spherically symmetric Fourier-transformed source function with simultaneous emission, 
Fij = fij{\Qij\) e ' ~" B ^*°, where the exponential term corresponds to emissions at a fixed time to. We believe that 
the assumption of instantaneous emission is a good approximation, because recent two-particle correlation data reveal 
a small time duration (i.e. R ut ~ -Rsidc)- In addition, the existence of a finite time duration would not affect the 
following treatment, because it does not influence the chaocitity in the two-particle correlation function but, rather, 
decreases the width of the correlation function. For the spatial part, fij, we consider a Gaussian, an exponential, and 
a function written in terms of cosh, which becomes quadratic at small |q| and exponential at large \q\: 

f 1 , ij (\q ij \)=exp(-R 2 \q ij \ 2 /2), (12) 
f 2 , ij (\q ij \)=exp(-R\q ij \/2), (13) 

/cosh(,R|qU) 



Here, q^ — p i — Pj, and R is the size parameter. The correlation functions are then calculated from the relations 

f 2 

C 2 (pi,p 2 ) = l + A inv -^f- (15) 
/11/22 



and 

\fjj) fain 711/22/33 / 

where Ai nv , v and v-$ are phcnomcnological adjustable parameters accounting for non-trivial coherence effects, and 
j) is a summation over (i, j) = (1, 2), (2, 3), (3, 1). The forms Eqs. (|15f) and (|16|) are chosen so as to represent the 
chaotic case for Aj nv = 1, v = 1 and 1/3 = 1. Equation (|15fl is a standard correlation function to fit data [l9(. Although 

an alternative form of the three-particle correlation function, C3 = 1 + v^P -\ -rH h 2v z l 2 { 12 { 23 { 31 , l is often used, 

HiSii J11/22/33 

this equation corresponds to a perfectly chaotic source with a long-lived resonance decay contribution, and is incapable 
of taking into account non-trivial coherence effects. Therefore, we simply parametrize the three-particle correlation 
function by including the parameters from both the two-particle correlation part [the second term in Eq. (|16|) ] and the 
genuine three-particle correlation part [the third term in Eq. I|16|l ] . Note that we need only C2 (0) — 1 and C3 (0) — 1 to 



1 Actually, the resultant \ 2 value for this expression is worse than that for Eq. 1161 . 
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obtain u>. We can set V3 = 1 without loss of generality for a description of C^Qz) at small Q3. Then, the parameters 
R, Aj nv and v are determined through a simultaneous x 2 fitting to both Ci{Q\2) and (^(Qa). 

The three-particle correlator r3{Q3)/2 is then obtained from separate projections of the numerator and denominator 
of Eq. (0} onto the single variable Q3 [2l| . The results are depicted in Figs. 0E] and listed in Table [I] 



TABLE I: Results of the \ 2 fitting to C 2 and C 3 



f(\q\) 


R [fm] 


A 


V 


X 2 /dof 




fi 


7.0±0.07 


0.54±0.01 


0.48±0.01 


110/30 


0.958±0.09 


h 


14.4±0.2 


1.18±0.03 


1.08±0.03 


79.7/30 


0.736±0.09 


h 


15.2±0.2 


0.71±0.01 


0.64±0.02 


15.8/30 


0.872±0.097 



The results show the failure of the Gaussian [Eq. (|12|) ] and exponential [Eq. 1|13[) ] fittings. In Figs [3] and |3J 
each curve looks in agreement with the data, but the values of x 2 /dof are larger than in the cosh source function 
case. Furthermore, the disagreement in ^(Qa) for the exponential (Gaussian) case reveals an over- estimation (under- 
estimation) of C2 in the extrapolation to small Qi nv . Constrastingly, the cosh source function [Eq. i|14|) ] yields excellent 
agreement with the experimental data. Note that the resultant errorbars of ±0.097 for u> in the cosh case can also 
be obtained by an extrapolation of the error band expressed by the dash-dotted curves in Fig. 0] The error band is 
calculated using the propagation law from the errors of R, Ai nv and v in the % 2 fit (1-er confidence level). Hence, our 
result for T3 in the cosh case is consistent with all the experimental data points, within the error bars. 

The deviation in Fig. 0] for larger Q\ is not significant. Because both the numerator and denominator in Eq. 
go to zero at the region, the projection requires a careful treatment and is strongly affected by unimportant tails. 
However, we do not know the real three-particle distribution which should be used as a weight. Therefore we use 
Eq. (|16|l instead. We confirmed that the most important quantity, u — ra(0)/2, does not depend on the choice of the 
weight. 



1 1 1 1 1 1 1 r 




10 20 30 40 50 60 70 80 90 
Q12 [MeV] 

FIG. 2: Two-pion correlation function C2{Q\2)- The curves represent the results of the fits for each source function. The dots 
represent the experimental results taken from Ref. [l9|. 

Now let us evaluate how chaotic the pion sources are by applying the models [Eqs. JSJ-JHl)] to the obtained values 
A truo = 0.93 ± 0.08 and ui — 0.872 ± 0.097, where we adopt the result of the cosh source function. The allowed 
parameter regions of Models I and II are summarized in Table [H] In these single-parameter models, e pc and a m 
should give consistent results for the two input quantities, A truc and u>. Table [n] reveals that, after the appropriate 
correction for the long-lived resonance decay contributions to the two-pion correlation chaoticity index, the two-pion 
correlation function gives results consistent with the three-pion correlation function. The value of e pc obtained from 
u> is slightly smaller than that obtained in the previous analysis carried out by the STAR collaboration 0] , while the 
values of e pc obtained from A true and ui are consistent with the result given in Ref. [l3| . The very large uncertainty on 
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FIG. 3: Three-pion correlation function C^Qz). Here, the curves and dots have the same identifications as in Fig. |5] The 
experimental data are taken from Ref. fl3|l . 
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FIG. 4: The Q3 dependence of r$/2. The curves and dots have the same identifications as in Figs. and |H] except for the 
dash-dotted curves, which represent the upper and lower bounds of the cosh case. The experimental data are taken from 
Ref. 0. See the text for details. 



a m is due to rapid changes of cv m as a function of A true and ui in these parameter regions. This feature may suggest 
the existence of a chaotic background introduced in Models I and III. 

The success of extracting the "true" chaoticity by considering the long-lived resonance decay contribution suggests 
that it is needed to determine the source of pions, i.e., direct and short-lived resonance decay or long-lived resonance 
dec ay. R ecently, a "partial" Coulomb correction has been proposed [23.1231 ] and actual computations have been carried 
out |24l l25j , in which pions from long-lived resonance decays are treated separately. Such a correction affects not only 
the HBT radii but also A exp . We leave its treatment for future publications |26j. 

The allowed region of e and a for Model III is displayed in Fig. In this figure, the darkest shaded area, labeled 
"III" , where the two lighter areas overlap, represents the ranges of e and a that correspond to the regions of Atruc 
and to estimated from the experimental data. The best fit point, e = 0.75 and a = 0.77, is indicated by the 
solid square. Unfortunately, the uncertainty is still too large to provide strong constraints on the pion production 
mechanism. Nevertheless, the result shows a strong correlation between a and e; that is, the experimental data reveal 
a mostly chaotic source which allows both a "partially coherent" picture (a ~ 0, s < 1) and a "multicoherent" picture 
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TABLE II: Results for Models I and II 



Model I Model II 



From A 


E pc = 0.73 ± 0.14 


Q m = 12.6 ± 14 


From uj 


e pc = 0.65 ±0.11 


a m = 7.6 ± 9.7 
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FIG. 5: Results for Model III. Area I is the allowed parameter region corresponding to the value of A indicated by the solid 
line (A truc = 0.85). Area II, which is bounded by the dashed curve (upper bound of ui) and dotted curve (lower bound of uj), 
is the allowed parameter region corresponding to the range of uj. Area III, which is characterized by overlap of areas I and II, 
is the allowed parameter region for a and e. The filled solid represents the best fit point. 



(a> l,e ~ 0). Though a fully chaotic source, e = 1, is not excluded, the strong Q 3 dependence of r 3 (Q 3 ) seen in 
Fig. m is inconsistent with chaotic and symmetric sources. Even if the source is perfectly chaotic, the asymmetry of 
the source could cause a strong Q3 dependence. In this case, however, r 3 becomes almost unity at small Q3, while 
it deviates from unity at large Q3 [Tlj • Hence, we can claim that small coherent components are more likely to be 
produced. 

In summary, we have given analyses of the two- and three-pion correlation data obtained by STAR, in terms of 
three kinds of partially coherent source models. We found that all three models give consistent description for both 
the three- and two-pion correlations if we take account of the apparent reduction of the chaoticity in the two-pion 
correlation function caused by the long-lived resonance decays. All models give the consistent picture that pion 
emission in the RHIC experiment is mostly chaotic, while coherent components are likely to exist (see Table II and 
Fig-EJ. Further interesting problems concerning the multiplicity and the collision energy dependence of the chaoticity 
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will be published elsewhere [2£ 
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APPENDIX A: DERIVATION OF A exp 



From Eq. J3J, assuming that the k t dependence of A cxp is governed by the long-lived resonance decay contributions, 
we can write down A^ xp for midrapidity as 

XT = mi " [ , I ■ (Al) 



where TVdir+short and Atotai represent the sum of the number of particles from direct emission and decay of short lived 
resonances and the total number of emitted particles, respectively. The values fcj . and fc( max are the minimum and 
maximum for a given momentum window of the z-th bin. Hence, k\ . = k\~~ . Then, the domain of integration for 
the momentum acceptance, k tmia = fc{ mIn < k t < fct' nax = &t max > with n being the number of the momentum bin, can 
be divided as 



fe 'mi„ 

In terms of the above relation, A P can be written as 



i=l 



(A2) 



Substituting Eq. IjAlj) into Eq. (|A3|). we obtain Eq. (TTJ 



E n i tfh* k t dk t ( dN t i '+t ott ) 
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